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ABSTRACT 

We revisit the bimodal distribution of the galaxy population commonly seen in the local 
universe. Here we address the bimodality observed in galaxy properties in terms of spec- 
tral synthesis products, such as mean stellar ages and stellar masses, derived from the ap- 
plication of this powerful method to a volume-limited sample, with magnitude limit cutoff 
M(r) = —20.5, containing about 50 thousand luminous galaxies from the SDSS Data Re- 
lease 2. In addition, galaxies are classified according to their emission line properties in three 
i/^ \ distinct spectral classes: star-forming galaxies, with young stellar populations; passive galax- 

■ ies, dominated by old stellar populations; and, hosts of active nuclei, which comprise a mix 
' of young and old stellar populations. We show that the extremes of the distribution of some 

^Q^i galaxy properties, essentially galaxy colours, 4000 A break index, and mean stellar ages, are 

I ■ associated to star-forming galaxies at one side, and passive galaxies at another. We find that 

O I the mean light- weighted stellar age of galaxies is the direct responsible for the bimodality seen 

^ . in the galaxy population. The stellar mass, in this view, has an additional role since most of 

^ ■ the star-forming galaxies present in the local universe are low-mass galaxies. Our results also 

^ ' give support to the existence of a 'downsizing' in galaxy formation, where massive galaxies 

seen nowadays have stellar populations formed at early times. 

■ Key words: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters 
^ ' — galaxies: stellar content — stars: formation 

■ 



1 INTRODUCTION 

It is usual to group galaxies in distinct types or classes, trying to 
find their general properties and to understand their differential for- 
mation and evolution. In the last few years, this galaxy 'taxonomy' 
has been interestingly simplified. As a by-product of recent redshift 
surveys, the study of galaxy populations has quantitatively revealed 
the existence of a bimodal distribution in some fundamental galaxy 
properties, found both in photometric and spectroscopic data. This 
bimodality of galaxy populations, despite its apparent simplicity, 
figures out as an important key in our whole comprehension of the 
processes which drive galaxy evolution. 

Perhaps the most representative bimodal distribution seen in 
galaxy properties is that found in galaxy colours. Since photometric 
parameters are easily measured, their bimodal behaviour has been 
studied for galaxies in the local universe by using both the Sloan 
Digital Sky Survey (SDSS) and the Two Degree Field Galaxy Red- 
shift Survey (2dFGRS) data iStrateva et all200ll : lHogg et all200l 
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iBlanton et all2003al : IWild et alJl2005h . and, for distant galaxies to 
z ^ 1, in the COMBO- 17 photometric redshi ft survev jBe ll et alJ 
l2004^ and in the DEEP galaxy redshift survey iWeiner et all2005h . 
This behaviour also persists for ev en higher redshifts (z ~ 1.4 
IWiegert. de Mello & Horelloil2004h . As galaxy colours reflect the 
star formation history of galaxies, a bimodality in their distributions 
suggests that they have evolved through two different major paths. 

Amongst local galaxies, iKauffmann et lll j2003d) have anal- 
ysed the star formation history and its dependence on stellar mass 
for a large sample of SDSS galaxies. They have found a transition 
in the physical properties of galaxies at a stellar mass ~ 3 x 
10^° Mq, which results in a bimodality between low-mass galax- 
ies with low concentration indices typical of discs and young stel- 
lar populations, and massive ones, with high concentration indices 
and old stellar populations. The connection between this bimodal 
distribution and tha t observed in galax y colours has been done, 
at least partially, by iBaldrv et"all i2004) in a work on the colour- 
magnitude diagram for a sample of SDSS galaxies. These authors 
have demonstrated that the transition in galaxy properties occurs 
around (1.5 - 2.2) x 10^° M© for the red sequence and around 
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(2 — 3) X 10^° M(T) for the blue seq uence, close to the value of 
found bv iKauffmann ^ <2003hh using spectroscopic data. This 
transition stellar mass is also associated to a shift in galaxy gas mass 
fractions, as pointed out by Kannappan (2004). Another impor- 
tant characteristic of bimodality in galaxy colo urs is its dependence 
on environment. This issue was addressed bv alogh etalJ i2004h 
who have complemented the analysis by Baldrv etaUiliOO^ to ac- 
count for, besides luminosity dependences, the environmental de- 
pendence of the red and blue galaxy fractions. 

A bimodal behaviour is also seen in the star formation prop- 
erties of galaxies. From the 2dFGRS data. lMadgwick et all J2002h 
have derived a spectral parameter 77 via a principal component anal- 
ysis method that is strongly correlated with the star formation rate 
(SFR) of a galaxy. A bimo dal distribution is clearly seen in this 
parameter JWild et alJl2005h . with a characteristic value that splits 
galaxies in two classes according to their star formation activity. 
Brinchmann et al. (2004) also find a similar distribution, showing 
that galaxies are from two kinds: one composed by concentrated 
galaxies, with low specific SFRs, and the other composed by less 
concentrated ones with high specific SFRs. 

In the first paper of this series on semi-empirical analysis of 
galaxies (Cid Fernandes et al. 2005a, hereafter SEAGal I) we have 
presented a spectral synthesis method able to derive the main phys- 
ical properties of galaxies from their spectral data only. The reli- 
ability of this approach was investigated by three different ways: 
simulations, comparison with other works, and an empirical anal- 
ysis of the consistency of results for a SDSS galaxy sample. We 
find that spectral synthesis provides reliable physical parameters, 
mainly stellar ages and stellar metallicities. Stellar masses, velocity 
dispersions and extinction are also obtained through this method. 
The analysis of residual spectra (observed minus synthesised) also 
provides useful information on the overall properties of emission 
line regions, like the nebular extinction and metallicities. 

The joint spectroscopic and photometric data provided by 
SDSS, in addition to the arsenal of physical parameters obtained 
by spectral synthesis, turns out as a great opportunity to study 
galaxy populations by analysing both spectral and structural clas- 
sifications. In this work, we will produce a spectral classification, 
initially based on the emission line properties of galaxies and the 
presence of nuclear activity, and subsequently we will compare 
this classification with a structural (morphological) one, which uses 
SDSS photometric information to separate galaxies between early 
and late-types. The bimodality of the galaxy populations will be 
analysed here with the starting point placed in these classifications. 

This paper is organised as follows. Section |2l presents an 
overview of our synthesis method, detailed in SEAGal I, focus- 
ing on the main outputs resulting from its application to a volume 
limited sample of SDSS galaxies. Section |3] describes the classifi- 
cation scheme used along this work, based on a diagnostic diagram 
of emission line ratios. Section |4| discusses the bimodal character 
of the galaxy population, with emphasis in the emission line prop- 
erties of galaxies. In Section|5]we discuss the main findings of this 
work. Finally, Section|6lsummarises our results. 



2 THE SPECTRAL SYNTHESIS METHOD 

In this section we present an overview of the spectral synthesis 
method used in this work and its main output. We also present the 
volume limited sample of SDSS galaxies analysed here, as well 
as our procedure for measuring emission line intensities from the 
residual spectra. 



2.1 The method 

In this work we use the STARLIGHT code to derive stellar popu- 
lation properties directly from the galaxy spectra (continuum and 
absorption lines only). STARLIGHT is built upon computational 
techniques originally developed for empirical population synthesis 
with additional ingredients from evolutionary synthesis models. In 
brief, we fit an observed spectrum with a combination of Simple 
St ellar Populations (SSP) fr om the evolutionary synthesis models 
of lBruzual & Cha riot ( 2003) (hereafter BC03). We use the "Padova 
1994" tracks recommended by BC03 and a Chabrier (2003) IMF 
between 0.1 and 100 M©. (Experiments with a Salpeter IMF with 
the same mass limits were also performed, yielding nearly identi- 
cal results except for stellar masses, which are 1.5 times larger than 
for a Chabrier IMF.) Extinctio n is modelled as due to foreground 
dust, with the reddening law of lCardelU. Clavton & MathiJ il98^ 
with Rv =3.1, and parameterised by the V-band extinction Ay- 
Line of sight stellar motions are accounted for using a Gaussian 
distribution. The fits are carried out with a simulated annealing 
plus Metropolis scheme, with regions around emission lines and 
bad pixels excluded from the analysis (see SEAGal I for details). 

Since SEAGal I, STARLIGHT has undergone a major revi- 
sion. The main technical change is that whereas in SEAGal I we 
worked with a single Markov chain all the way from an initial 
guess to the final model, we now run several parallel chains. At 
each level of the annealing scheme (i.e., for each "temperature") 
we run the chains u ntil they satisfy a converg ence criterion similar 
to that proposed bv iGelman & Rubin il992h . Adaptive step sizes 
are employed to improve efficiency. These and other details will be 
described in a future communication accompanying a public ver- 
sion of STARLIGHT. 

A second and more important difference with respect to SEA- 
Gal I is that we now use a larger base of SSPs, comprising — 
150 elements of 25 different ages between 1 Myr and 18 Gyr, and 
6 metallicities: Z = 0.005, 0.02, 0.2, 0.4, 1 and 2.5 Zq. The age 
grid was chosen by exploring an implicit assumption in spectral 
synthesis with a discrete base, namely, that the spectrum of a pop- 
ulation with age t between tj and t^+i is well represented by a 
linear combination of the j and j + 1 SSPs. We started from the 
set of 15 ages in the SEAGal I base, and then fitted each of 221 
BC03 SSP spectra for the same Z with the two components which 
bracket its age. This allows us to identify SSPs whose spectra are 
poorly interpolated by combinations of the closest elements in the 
base, and thus should be included in a finer base. Following this 
procedure (graphically illustrated in Leao 2006), we arrived at a set 
of 25 ages. The difference between original and interpolated SSP 
spectra with this base is better than 1 per cent on average, exceed- 
ing 5 per cent only for a few t < 10^ yr models at Z ^ O.2Z0, 
corresponding to short-lived phases when massive red-supergiants 
show up. 

The refinement in the age grid is not as relevant as the in- 
clusion of lower metallicity components in the base. Whereas the 
— 45 base in SEAGal I comprised only Z ^ 0.2Zq, we now 
allow for Z = 0.005 and 0.02^©. These very metal poor SSPs 
are bluer than those of Z ^ 0.2Zq SSPs, particularly at old ages. 
Such low Z populations, if present in a galaxy but ignored in the 
fits, may lead to an underestimation of ages, as the code tries to 
compensate for their lack with younger components of higher Z 
due to the age-Z degeneracy. To illustrate this point. Fig [I] shows 
the predicted evolution of the 4000 A break index 1)^(4000) (Sec- 
tion]^) and its observed distribution for the sample described be- 
low. The peak at i^n(4000) ^1.3 matches well the value reached 
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Figure 1. (a) Observed distribution of the 4000 A break index for our vol- 
ume limited sample, (b) Evolution of Dn(4000) for SSPs of 6 different 
metallicities. From right to left: Z = 2.5, 1,0.4, 0.2, 0.02 and 0.005 Zq.(c) 
Relation between their Dn(4000) and z-band mass-to-light ratio for SSPs 
of 6 different metallicities. From right to left: Z = 2.5, 1, 0.4, 0.2, 0.02 
and 0.005 Zq . The dotted lines mark the location of the lower peak in the 
observed Dn(4000) distribution. 



by Z ^ O.O2Z0 SSPs of ages from 1 Gyr to over 10 Gyr. To re- 
produce this value with Z ^ 0.2^© one must either accept that 
these galaxies formed less than ^ 1 Gyr ago or else invoke more 
continuous star-formation regimes (exponentially decaying, multi- 
ple bursts or other variants). Whereas the latter is in fact a likely 
explanation for galaxies with low Dn(4000), our point here is that 
by excluding low Z components from the base one would be in 
practise imposing this interpretation a priori. Besides this "math- 
ematical" argument, the mere fact that very metal poor stars exist 
provides physical motivation to include them in the base. Having 
said that, one must also recall that the synthetic SEDs of such pop- 
ulations are affected by a number of caveats, as discussed by BC03. 

One must also realise that the inclusion of 0.02 and 0.005 Zq 
SSPs is bound to have an effect on the stellar masses. This is illus- 
trated in the top panel of Fig[l] which shows that the mass-to-light 
ratio in the z-band of ^ 10 Gyr, 0.02 and 0.005 Zq SSPs is up to 
one order of magnitude larger than for Z ^ O.2Z0 models of simi- 
lar Dn (4000). To monitor the magnitude of this difference we have 
also fit the data with a base excluding the 0.02 and 0.005 Zq SSPs. 
Throughout the paper we will mention results obtained with this 
alternative base whenever they differ substantially from the ones 
obtained with the full set of metallicities. We anticipate that such 
changes are mostly quantitative; they do not affect the overall qual- 
itative conclusions of our analysis. 

Another difference with respect to the base in SEAGal I is that 
we now allow for SSPs as old as 18 Gyr, which is formally incon- 
sistent with the ^ 13.5 Gyr old Universe implied by our adopted 
cosmology, ifo = 70 km s~^ Mpc~\ Qm — 0.3 and Qa — 0.7. 
We relax the cosmological limit merely to avoid the crowding of 
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Table 1. Summary of parameter uncertainties. Each row corresponds to a 
parameter constructed by combining the parameters given by the synthe- 
sis. The different columns list the rms difference between output and input 
values of the corresponding quantity, as obtained from simulations with dif- 
ferent signal-to-noise ratios (see SEAGal I for details). The units are dex for 
logarithmic quantities and mag for Ay . 

mean stellar ages at the upper age limit (visible, for instance, in 
SEAGal I's Fig. 13b). Given the uncertainties in stellar evolution, 
cosmology, and the spectrophotometry measurements, we do not 
regard this inconsistency as critical, and, in any case, we stress that 
this choice has no consequence for the main conclusion of this pa- 
per. 

The result of this procedure is a list of parameters for each 
galaxy. As explained in SEAGal I, instead of working with the 
full population vector, whose individual components are hopelessly 
plagued by mathematical and astrophysical degeneracies, we de- 
scribe the results of the spectral synthesis in terms of a small num- 
ber of robust parameters. The most important for the discussion 
in this paper are: a) the dust-corrected present stellar mass inside 
the fibre; the total stellar mass (M^) is computed a posteriori after 
correcting for the fraction of luminosity outside the fibre assuming 
that the galaxy mass-to-light ratio does not depend on the radius 
(but see Appendix A); b) the mean stellar age, weighted either by 
light, (logt^)L, or by mass (log t^)M; c) the mean stellar metal- 
licity, light and mass- weighted, log(Z^) l and log(Z^)M, respec- 
tively; and d) the V-band stellar extinction Ay. A summary of the 
uncertainties in the main parameters that will be used in this work 
is shown in Table [I] These uncertainties were determined in SEA- 
Gal I through the difference between output and input values of the 
corresponding quantity, as obtained from simulations with different 
signal-to-noise ratios. In addition to this list of parameters, by sub- 
tracting the synthesised spectrum from the observed one, we get a 
"pure-emission" residual spectrum, which is useful for analysis of 
the galaxy emission lines. These are discussed next. 

2.2 Application: SDSS data 

We apply our synthesis method to a large sample of SDSS galaxies 
aiming to investigate the properties and the origin of the bimodal 
galaxy distribution. The SDSS spectra cover a wavelength range 
of 3800-9200 A, have a mean spectral resolution A/ A A 1800, 
and were taken with 3 arcsec diameter fibres. The spectroscopic 
sample studied here is the same discussed in SEAGal I. Briefly, we 
select ed a volume liniited s ample from the SDSS Data Release 2 
rDR2 lAbazaiian et alJl2004 . with a redshift range of 0.05 < z < 
0.1, corresponding to an absolute magnitude limit cutoff M(r) = 
—20.5, or M*(r) + 1. The absolu te magnitudes adopted here are 
/c-corrected iBlanton et alJl2003bl) . We also restricted our sample 
to objects for which the observed spectra show a signal-to-noise 
ratio greater than 5 in the ^, r, and i bands. In addition, we have 
detected a minor fraction of galaxies with multiple spectral data. 
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from which we have excluded those with smaller signal-to-noise 
ratios. All these criteria leave us with a sample containing 49917 
galaxies, corresponding to a completeness level of 98.5 per cent. 

The spectra are corrected for Galactic e xtinction with the maps 
of bv ISchleg eL Finkb einer & DavisI il998h and the extinction law 
of ICardelli et al.l ([l283), shifted to the rest-frame and resampled 
from 3400 to 8900 A in steps of 1 A. Notice that this is a sHghtly 
wider range than that employed in SEAGal I (3650-8000 A), al- 
though not all galaxies have data close the edges of this new spec- 
tral range. 

2.2.7 Masks 

A further novelty is that now each galaxy has its own emission line 
mask, as opposed to the single general mask for all galaxies used 
in SEAGal I. General masks may unduly throw away data when 
emission lines do not exist, and, conversely, ignore lines which are 
normally weak but may occasionally be strong enough to deserve 
being masked, (e.g. He II A4686). The individual masks are con- 
structed from an initial STARLIGHT fit with a general mask. The 
strongest lines are fitted in the residual spectrum (observed minus 
model) using the code described in the next section, and the infor- 
mation on line widths and velocity off- sets are used to search for 
other emission lines out of a list of 52 lines in the 3400 to 8900 
A range. When present, these lines are masked in windows whose 
sizes are defined by the A's at which their assumed Gaussian pro- 
files become weaker than 1/2 the local rms noise flux density. Vi- 
sual inspection of hundreds of cases shows that this scheme is able 
to identify and mask even weak emission lines satisfactorily. 

Besides emission line masks, we also exclude four other win- 
dows from the fits: 5880-5906 A, to skip the Na D AA5890,5896 
doublet, which is partly produced in the interstellar medium; the 
6850-6950 and 7550-7725 A, for which BC03 had to resort to the- 
oretical spectra due to problems in these ranges in the STELIB li- 
brary. Curiously, when residual spectra are averaged over thousands 
of galaxies, these two windows st and out as stro ng "emission" fea- 
tures which are clearly spurious fGomes"2005^. The last window, 
in the 7165-7210 A range, shows a similar systematic broad resid- 
ual in emission, so we have masked it out too, even though it is not 
listed as problematic by BC03. 



2.3 Emission line measurements 

One of the products of our spectral synthesis procedure is a resid- 
ual spectrum suitable for measurements of emission lines, as the 
stellar continuum and absorption features are well modelled. Even 
weak lines, as the [O lll]A4363 auroral line, can be easily mea- 
sured (when present, of course ). Indeed, such an app roach has been 
adopted by other authors (e.g. iTremonti et alJbOoll) to investigate 
SDSS emission line spectra. 

We have developed a code to measure the intensity of main 
emission lines from the residual spectrum by fitting Gaussian 
functions to the line profiles, characterised by three parameters: 
width, offset (with respect to the rest-frame central wavelength), 
and flux. Lines from the same ion are assumed to have the same 
width and offset. Additionally, [O lll]A5007/[O lll]A4959 = 2.97 
and [N ll]A6584/[N ll]A6548 = 3 flux ratio constraints are im- 
posed. Furthermore, only lines detected with signal-to-noise ra- 
tio (S/N) greater than 3 are conside red in the an alysis. This is 
the same value adopted bviBrinchman n et alli2004h in their study 
of physical properties of star-forming SDSS galaxies, and is la 
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Figure 2. The BPT diagram for our subsample of galaxies with emission 
lin es measured with S/N gr eater than 3. The solid line is the curve defined 
b^jKauffinam^^^ J2003cl) and the dashed one is the curve defined by 
lKew1evetalJi200ll) (see text for details). In this and other figures, the lin- 
ear grey scale level represents the number of galaxies in each pixel, darker 
pixels having more galaxies. 

greater than the detection limit adopted by iMiller et all i2003h 
in an environmental analysis of active galactic nuclei, also in 
the SDSS. The lines that are currently being measured by our 
code include [O ll]AA3726,3729, H^, H7, [O lll]A4363, H/3, 
[O iii]AA4959,5007, [O i]A6300, [N ii]A6548, Ua, [N ii]A6584 
and [S II] AA67 16,6731. For each emission line, our code returns the 
rest-frame flux and its associated equivalent width (EW), the line 
width, the velocity displacement relative to the rest-frame wave- 
length, and the 5'/A^ of the fit. In the case of B aimer lines, the un- 
derlying stellar absorption is also measured directly from the syn- 
thetic spectra, yielding absorption fluxes and equivalent widths. 

Additionally, we have also identified broad emission 
lines in the residual spectra, characteristic of galaxies host- 
ing Seyfert 1 nuclei, following a proc edure similar to that of 
IVeron-Cettv. Veron & GoncalveJ i200lh . As a result of this anal- 
ysis we have identified 335 objects with significant broad compo- 
nents. 



3 DEFINITION OF SPECTRAL CLASSES OF GALAXIES 

In this section we present the galaxy types that will be used here- 
after. This galaxy classification is based on spectral properties, i.e. 
on the presence or absence of emission lines and on their ratios in 
the case they are present in the spectra. 

3.1 Galaxy classification 

In general, the classification of narrow emission line galaxies is 
made according to the mechanism by which lines are produced, 
grouping galaxies into two types. The first group is composed 
by normal star-forming galaxies, with H II region-like spectra, in 
which the gas is photoionised by young, hot OB stars. The second 
one is formed by galaxies hosting active galactic nuclei (AGN), 
which produce a much harder radiation field ionising the surround- 
ing gas. 

These emission line galaxies are generally classified with 
the help of diagnostic diagrams formed by line ratios of the 
strongest lines present in their spectra. Here we will use the 
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Figure 3. Examples of galaxy spectra of the spectral classes analysed in this 
work. The residual spectra are also shown as thin lines at the bottom of each 
panel. These spectra have been smoothed in order to improve appearance. 



[O III1A5 007/Hg versus [N ll1A6584/Ha dia gram firstly pro- 
posed by Baldwin. Phillips & T erlevicl] il98lh and updated by 
IVeilleux & Osterbrocb Il987i) : hereafter we will refer to it as the 
BPT diagram. In it, galaxies form two distinct branches, which look 
like the wings of a seagull. The left wing contains star-forming 
galaxies and the observed sequence corresponds to a change in 
metallicity of the H II regions emitting the lines. The right wing 
appe ared clearly only with the most recent galaxy surveys (e.g. 
iKau ffmann et al. 2003c) and corresponds to galaxies hosting ac- 
tive nuclei; photoionization models show that galaxies in this wing 
cannot be ionised only by radiation from massive stars, an addi- 
tional heating/i onising source is ne cessary to explain the observed 
line ratios fe.g. lKewlev et all200lh . 

The BPT diagram for our sample is shown in Fig.|2| We have 
considered only those galaxies having a signal-to-noise ratio in all 
four emission lines greater than 3; this limit results in a subsample 
containing 23080 narrow emission line galaxies. Here we have not 
considered galaxies with spectra containing broad emission lines 
(following Secti on_ l23l a bove). After an analysis of SDSS data, 
IKauffmann et alJ (20033) have defined a curve (shown in Fig.lJJas 
a solid line) in the BPT diagram to distinguish pure star-forming 
galaxies from those contaminated by nuclear non thermal emission. 
Although the exact position of this curve could be questionned (see 
^Brinchman n et al . 2004 ; Stasins ka et al. 200^, we hav e adopted 
it for easy reference to recen tly published work (e.g. iHao et all 
I2OO5I : Ifirinchmann et al JbOoH) to distinguish normal star-forming 
galaxies (below the curve) from galaxies hosting an active nucleus 
(above it). In Fig.|2|we also show the curve given bv lKewlev et alJ 



i2001l) . which we use in Sect. 4.2.1 to define galaxies with 'com- 
posite' spectra, as done bv lBrinchmann et a'll i2004). 

The normal star-forming galaxies in our sample show sig- 
nificant amount of star formation (more than 99 per cent have 
EW(Ha) > 5A). Note that as we are studying a volume limited 
sample built with an absolute magnitude cutoff M(r) < —20.5, in 
this diagram we do not see the faint star-forming galaxy population 
that inhabits the extremes of its left wing (corresponding to a line 
ratio of [N Il]/Ha < 0.1). 

Galaxies hosting active nuclei (hereafter AGN hosts), on the 
other hand, populate the upper right wing of the BPT diagram. In 
addition, we do not distinguish these galaxies from the 'composite' 
or 'transition' galaxies, which have a substantial contribution of 
their Ha emission due to active nuclei, and populate the region 
near the solid curve shown in Fig. 121 

The objects that do not appear in Fig. |2l because of insuffi- 
cient S/N in at least one of the four relevant lines also receive a 
classification. For instance, galaxies that host active nuclei are also 
identified if they have only the e mission lines of Ha and [N ll] 
measured with 3cr confidence fe.g. lCoziol et alJll998l:lMiller et alJ 
l2003h . We have selected these galaxies by considering the limit 
log ([N ii]/Hq;) > -0.2, and classified them as AGN hosts. We 
note that the fraction of galaxies with active nuclei identified by 
this way is about 54 per cent of the total of AGNs hosts in our 
sample. Additionally, galaxies with broad emission lines which 
host Seyfert 1 nuclei, are naturally added to the spectral class of 
AGN hosts. Galaxies without evidence of significant star forma- 
tion, whose spectra do not show the emission lines of Ha and H/3, 
or show them with equivalent widths smaller than 1 A, are classi- 
fied as passive galaxies. Note that with this definition we are not 
including the fraction of elliptical galaxies which show eviden ce of 
star f ormation activity fe.g. iNakamura et al. 2004 : Fukugita et alJ 
l2004h in our subsample of passive galaxies. The remaining galax- 
ies, for which we can not give a classification based on the BPT 
diagram, tend to include mainly star-forming and passive galaxies. 
These unclassified galaxies represent a reduced fraction of galaxies 
in our sample (about 3 per cent) and will not be analysed in this 
work. 

A summary of the spectral classes described above, which will 
be analysed in the following sections, is shown in Table |2| The 
number and per percentage of each spectral class, as well as the me- 
dian values of their concentration indices, colours, 4000 A break, 
mean light-weighted stellar ages and stellar masses, are shown in 
this table. In Fig.|3]we show some examples of galaxy spectra for 
the distinct classes of galaxies discussed above. In each panel is 
also shown the residual spectrum after removing the model ob- 
tained by our synthesis method. These spectra correspond to galax- 
ies with both stellar age and stellar mass equal to their median val- 
ues for each class, therefore they are very illustrative of the different 
spectral features which characterise the classes we discuss here. 

3.2 Sources of bias 

As mentioned in Section l2!3l in this work we are imposing a limit 
in the S/N of the measured lines in order to consider them in 
emission or not. The choice of this limit affects the distribution 
of objects in each spectral class, as we shown in Table |3| where 
the number and percentage of galaxies in our classes are listed for 
three limiting values of the S/N. For a change in S/N going from 
2 to 5, the proportion of star-forming galaxies increases by only 6 
per cent, but that of passive galaxies increases by 50 per cent. In 
contrast, the proportion of active nuclei hosts decreases by about 
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Class 


Number 


Per cent 


C 


iu-i) 






log M^/M0 


All 


49917 


100.00 


2.67 


2.76 


1.69 


9.53 


10.79 


Star-forming 


16108 


32.27 


2.25 


2.10 


1.35 


8.91 


10.56 


Passive 


10485 


21.00 


2.97 


3.06 


1.93 


9.86 


10.87 


AGN Hosts 


21733 


43.54 


2.77 


2.88 


1.78 


9.66 


10.90 


Unclass. 


1591 


3.19 


2.67 


2.83 


1.73 


9.58 


10.77 



Table 2. Summary of the spectral classes defined in this work and the median values of their mean light-weighted stellar ages and stellar masses. 



S/N 


Star-forming 


Passive 


AGN Hosts 


Unclass. 


2 


16024 32.1 


8626 17.3 


24321 48.7 


946 1.9 


3 


16108 32.3 


10485 21.0 


21733 43.5 


1591 3.2 


5 


17034 34.1 


12998 26.0 


17382 34.8 


2503 5.0 



Table 3. Number and percentage of objects in each spectral class as a func- 
tion of the limit in S/N when considering a line in emission. Here we are 
using 2L S/N limit equals to 3. 



S/N limits 


Star-forming 


Passive 


AGN Hosts 


Unclass. 


2 and 3 


99.99 


99.99 


50.97 


31.38 


3 and 5 


99.92 


99.99 


0.30 


25.12 



Table 4. Kolmogorov-Smirnov (KS) probability parameter for redshift dis- 
tributions when comparing the S/N limits of 2, 3 and 5 for considering a 
line in emission. Results are shown for each spectral class. Low values of 
the KS probability imply statistically different distributions. 

28 per cent. It is interesting to note that the fraction of unclassified 
galaxies is always smaller than 5 per cent. 

We have investigated how the choice of the S/N limit affects 
our analysis by looking for differences in the redshift distributions 
of our classes. This approach can also test possible aperture bias 
in our class definition. Since SDSS spectra are taken with only 3 

— , — , — , — , — I — < — < — < — , — I — , — , — , — , — I — , — , — , — I — I — I — I — , — ^ 

star-forming 

0-6 - Passive 

AGN hosts 




0.05 0.06 0.07 0.08 0.09 0.1 

redshift 

Figure 4. Testing aperture effects. The fraction of galaxies in each spec- 
tral class is shown in different bins of redshift. Solid lines are for fractions 
obtained with the S/N > 3 limit for detection of emission lines, whereas 
dashed lines are for the limit S/N > 5. The error bars follow a Poisson 
statistics. 



arcsec fibre diameter, for the nearest galaxies they might tend to 
sample mainly the bulge component, increasing the fraction of non 
emission line galaxies (especially passive ones). Therefore, the red- 
shift distributions for each spectral class defined by using different 
S/N limits, will be similar if the aperture bias does not play a sig- 
nificant role in our sample, or it does in an uniform way indepen- 
dently of the S/N value. In Table |4] we show the Kolmogorov- 
Smirnov (KS) probability parameter for the redshift distributions 
of our distinct spectral classes. We have compared the distributions 
of classes defined by using a S/N limit equals to 2 and 5, with 
those of our classes defined by the adopted value of S/N = 3. 
The values of the KS probability are very large for passive and 
star-forming galaxies, indicating statistically similar distributions 
for these spectral classes. In the case of AGN hosts, the value of 
the KS probability parameter is large when comparing S/N limits 
of 2 and 3, but it is very small for 5'/A^ limits of 3 and 5, indicating 
statistically different redshift distributions when comparing these 
S/N limits. For unclassified galaxies the KS probability values are 
larger than 25 per cent in the two cases investigated here. 

In Fig. 131 we show the redshift dependence of the fraction of 
galaxies in each spectral class obtained by considering two values 
of the S/N limit to detect emission lines (3 and 5). The expected 
effects due to aperture bias would tend to, for instance, increase the 
fraction of star-forming galaxies relative to that of passive along 
the redshift range of our study. In fact, what we note in this figure 
is an increment in the fraction of star-forming galaxies with in- 
creasing redshift. In addition, for lower redshifts, we also observe a 
higher fraction of AGN hosts probably as a result of the significant 
sampling of the central regions of nearby galaxies. In contrast, the 
fraction of passive galaxies does not significantly increase in the 
regime of low redshifts. Following these trends, we infer that our 
classes seem to be affected by aperture bias in the sense that the 
fraction of star-forming galaxies increases with redshift, and that 
of AGN hosts decreases with z. Some aspects of this effect in our 
sample are discussed in Appendix IXI 



4 BIMODALITY OF THE GALAXY POPULATION 

In this section, we discuss the characterisation of the two galaxy 
populations that inhabit the local universe. Our starting point is the 
definition of the spectral classes discussed in the last section. 

4.1 Spectral versus structural galaxy types 

Many works on the SDSS data have found that the concentra- 
tion index (C), defined as the ratio of Petrosian Rqo to R50 
radi i in r-band can be used to separate early and late-type galax - 
ies JShimasaku etal]l200ll : lstrateva et alJboQllilGoto et alJl2002h . 
since early-type galaxies tend to have light profi l es more centrally 
concentrated than the late-types iMorganl[l958h . lShimasaku et alJ 
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(c) Passive 
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Figure 5. Concentration index versus {u — r) colour for all galaxies in our sample (panel a) and for our distinct spectral classes (panels b-d). In this and 
in other figures of this paper, the top numbers on the right are the number of galaxies in each panel and the Spearman rank correlation coefficient (r^). The 
median relation between colour and concentration index in three bins of galaxy r-band luminosity containing the same number of galaxies is shown as distinct 
lines. The ranges of each bin are shown in the top-left legend. 



i200lh studied this parameter for a sample of morphologically clas- 
sified bright galaxies from SDSS. They found a strong correla- 
tion between C and morphological types, suggesting that this pa- 
rameter can be used to classify galaxies morphological ly (see also 
|Doi, Fukugita, & Okamura|[T993|; I Abraham et aljlB- Howeyer, 
they also noted that it is difficult to construct a pure early-type 
galaxy sample based only on the concentration index, since the 
resulting sample has ^ 20 p er cent contamination by late-type 
galaxies. IStrateya et alJ i200lh also studied the reliability of the 
concentration index to separate early and late-type galaxies. They 
adopt a yalue of C = 2.63 as a morphological separator (differ- 
ent from Shimasaku's yalue of C ^ 3). On the other hand, since 
galax y colours are a more conyentional estimator of galaxy types 
(e.g. JSandagd Strateya et al. concentrate their analysis on 

the (u — r) colour, finding that (u — r) = 2.22 clearly separates 
early (E, SO, and Sa) and late (Sb, Sc, and Irr) morphological types, 
with the {u — r) colour also expected to correlate with Hubble type. 

How are these parameters associated with the spectral classes 
we defined before? In Fig.|5]we show the relation of these two main 
morphological identifiers, C and (u — r), for all objects in our sam- 
ple and discriminated according to our distinct spectral classes (ex- 
cept for unclassified galaxies). In this figure, we also plot as lines 
with different shapes the median yalues of {u — r) colour along 
bins of concentration index for three luminosity ranges contain- 
ing the same numb er of galaxies. As found by other works (e.g. 
IStrateya et alJbOOlT) . the concentration index and {u — r) colour 
present an eyident correlation (quantified by the Spearman non- 
parametric correlation coefficient, rs — 0.65) when all objects are 
considered. Star-forming galaxies also show a correlation between 
these parameters, in a less significant leyel, in the sense that red 
star-forming objects tend to be more concentrated than blue ones. 
On the other hand, passiye galaxies do not show a correlation; the 
colours of these objects are almost independent of their concentra- 
tion, in all ranges of galaxy luminosity. Indeed, for passiye galaxies 
in our sample we haye the following mean yalues and dispersions: 
(C) = 2.96 =b 0.08 and {u - r) ^ 2.71 ± 0.22. In Fig.glwe also 
note that star-forming and passiye galaxies are responsible for the 
distinction of the two main groups of galaxies in the colour-C dia- 
gram. The class of AGN hosts is composed by mixtures of blue and 
red galaxies, with a large range in the concentration parameter. 

We show in Fig.|6t-b the distributions of concentration index 
and (u — r) colour for our spectral classes (also including unclas- 
sified galaxies). In these figures, we can easily see the dominance 
of star-forming galaxies in the blue and less concentrated galactic 
population; passiye galaxies, contrarily comprise the reddest and 



more concentrated galaxies. We also note that AGN hosts occupy 
an intermediate locus in both colour and concentration index distri- 
butions, showing again the mix of populations present in this spec- 
tral class. 

From the distributions shown in Fig.|6l we can obtain an opti- 
mal yalue that separates those two extreme classes formed by star- 
forming and passiye galaxies. Following . Strateya et al. (2001). we 
can define two parameters, reliability and completeness, which are 
used to do this task. The reliability (T^sf and T^p, for star-forming 
and passiye spectral classes, respectiyely) is the fraction of galax- 
ies from a giyen spectral class that are correctly classified by us- 
ing the optimal yalue, while the completeness (Csf and Cp) is the 
fraction of all galaxies of a giyen spectral class that are actually 
selected. In this way, we can find the optimal yalue that maximise 
the product CsfIZsfCpIZp (c.f. Baldry et al. 2004). We find that 
C — 2.62 and {u — r) — 2.35 are the optimal separators among 
star-forming and passiye galaxies (close to the yalues obtained by 
Strateya et al. to classify galaxies into early and late-types), with re- 
liability and completeness parameters (in per cent) of T^sf = 86.0, 
7^p = 89.8, CsF = 93.0 and Cp = 80.3 for concentration in- 
dex and 7^sF = 93.6, 7^p = 93.8, Csf = 96.0 and Cp = 90.4 
for {u — r) colour. Here it is important to stress that we are not 
accounting for the yariation of the colour diyi der along the magni- 
tude range of our sample. Indeed, as shown by lBaldry et al.N2004h 
the colour yalue which diyides the blue and red populations should 
slightly decrease with decreasing galaxy luminosity. 

Another parameter useful to distinguish between early and 
late-type galaxies is the 4000 A break, which is small for galax- 
ies with younger stellar populations, and large for older galaxies 
(as shown in Fig. 13 of SEAGal I). We me asured this index fol - 
lowing the narrow definition introduced by iBalogh et al.l il999h . 
with the continuum bands at 3850-3950 and 4000-4100 A; this 
narrow index will be referred as Dn(4000). Fig.|6j: shows the dis- 
tribution of this index for our spectral classes. We clearly note 
the bimodal distribution of this parameter, with a diyisory line at 
i^n(4000) = 1.67 separating star-forming and passiye galaxies 
with the highest yalues of both reliability and completeness (> 98 
per cent). 

In Fig. we show the Dn(4000) index as a function of both 
concentration index and (u — r) colour for all galaxies in our sam- 
ple, and for the spectral classes analysed here. The median yal- 
ues of i^n(4000) in bins of C and {u — r) for three ranges of 
galaxy luminosity are also shown as lines with different shapes. 
When all galaxies are considered, the Dn(4000) index shows a 
correlation with both concentration index {rs = 0.65) and colour 
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Figure 6. Distributions of (a) concentration index C, (b) {u — r) colour, and (c) Dn (4000) index, for all galaxy sample, star-forming, passive and AGN hosts 
spectral classes; the distributions for unclassified objects are also shown. The vertical dashed lines are the best population separators when we consider only 
distributions of star-forming and passive galaxies. 



{rs — 0.81). For star-forming galaxies there is a significant corre- 
lation only with colour {rs — 0.63), implying that these two quan- 
tities are linked for this type of galaxies. A natural explanation for 
these relations is that young and blue stellar populations are respon- 
sible for maintaining the correlation among colour and Dn (4000) 
for star-forming galaxies. On the other hand, as the concentration 
index is more associated with the shape of a galaxy than with its 
stellar population properties, the absence of correlation among C 
and Dn (4000) index reflects that star-forming galaxies with same 
mean stellar age (or Dn (4000) value) span a wide range of mor- 
phological types (or concentration values). We also note in Fig.0 
that passive galaxies do not show significant correlations among 
the quantities analysed here, and galaxies hosting AGN show cor- 
relations in both colour and C versus Dn (4000) plots, as expected 
due to the mixture of populations comprised by this spectral class. 

We thus conclude that by using spectroscopic information 
(emission lines and spectral features) one can obtain a better split- 
ting between the star-forming and passive galaxy populations than 
by using concentration and colours. Our results also indicate that 
galaxies hosting AGNs tend to present intermediate behaviour with 
respect to these two populations. Additionally, hereafter we will 
consider the star-forming and passive galaxies as the extremes 
of both blue and red galaxy distributions. Moreover, the optimal 
values that distinguish these spectral classes, mainly that for the 
Dn (4000) index, will be used to define the two main galaxy pop- 
ulations which inhabit the local universe, historically the early and 
late galaxy types. 



4.2 Physical properties of spectral classes 

Here we analyse the properties of the main physical parameters 
derived from our spectral synthesis method for each spectral class 
defined above. 



4. 2. 1 Mean stellar age 

The mean light- weighted stellar age of a galaxy reflects the epoch 
of formation of massive and bright O and B stars, frequently asso- 
ciated to starbursts. Thus, it is strongly affected by the recent star 
formation history of a given galaxy. On the other hand, the mean 
mass-weighted stellar age is associated with the epoch of formation 



of the stellar population which nowadays contributes most signif- 
icantly to the galaxy mass. Consequently, it is associated with the 
mass assembly history of a given galaxy. 

In Fig. [S] we show the distributions of these two age estimates 
for all galaxies in our sample, and for each spectral class separately. 
The distribution for the mean light-weighted stellar age shows a 
conspicuous bimodality with the star-forming and passive galaxies 
at their extremes, and the spectral class of AGN hosts occupying an 
intermediate locus in it. For this distribution, (log t^)L 9.53 is 
the age which best divides the extreme galaxy populations (shown 
in Fig.[8]as a vertical dashed line). On the other hand, the distribu- 
tion of (log t*)M is single peaked at older ages, with median value 
of about 12.4 Gyr for all sample. This implies that galaxies of all 
spectral classes in our sample have a large fraction of their stellar 
masses in old stellar populations. Thus, it is seen that the bimodal- 
ity noted in some galaxy properties, as discussed above, is pri- 
marily related to 'light- weighted' quantities. In other words, recent 
episodes of star formation present in star-forming galaxies, in con- 
trast with the absence of significant amount of young stars in pas- 
sive galaxies, produce the bimodal distribution seen in the (log t^) l 
distribution, as well as in the {u — r) colour and Dn (4000) index 
distributions. 

We have also investigated the behaviour of these trends for the 
results obtained with a base excluding the low-Z SSPs, which, as 
shown in Fig.[l] tend to result in older ages for young stellar popu- 
lations. Actually, we found this is true: the values of the mean stel- 
lar ages obtained by using the base without the 0.005 and 0.02 Zq 
SSPs, compared to that obtained with the base adopted in this work, 
decrease about 0.18 dex when all galaxies are considered and 0.31 
dex for star-forming galaxies. The age with best divide the galaxy 
populations for this 'metallic' base is (logt^)L — 9.32. We also 
note that the distribution of (log t^)M for this base shows a second 
peak at younger ages, besides that at older ages. As discussed in 
Section im this is a consequence of using Z ^ 0.2^© in the spec- 
tral base, which imply that one has to assume that younger galaxies 
were formed less than ^ 1 Gyr ago, contrarily to the assumption of 
more continuous star-formation regimes generally invoked to ex- 
plain the properties of these galaxies. 

As shown in Fig.jsj AGN hosts constitute an intermediate pop- 
ulation of galaxies. We analyse this behaviour of AGN hosts by in- 
vestigating the distribution of the [O lll]A5007 lu minosity, L\0 III] 
(thoug ht to be a tracer of the AGN power, e.g. iKauffmann et alJ 
l2003ch . as a function of the mean light-weighted stellar age. This 
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Figure 7. Concentration index and (u — r) colour versus Dn(4000) for all galaxies in our sample and for our distinct spectral classes. It is also shown the 
number of galaxies in each panel as well as the Spearman rank correlation coefficient (r^-). The median values of Dn(4000) in bins of C and {u — r) for 
three ranges of galaxy luminosity are also shown as lines with different shapes. 




relation is shown in Fig.|9l where we also plot the median values of 
L[0 III] as a function of (log t^) l for three ranges in galaxy lumi- 
nosity (shown as lines with different shapes); the age value which 
best separates star-forming and passive galaxies is also shown as 
a vertical dashed line. The correlation among these quantities is 
in the sense that the [O III] luminosity is low in older galaxies, 
with median values larger for brighter galaxies. In galaxies hosting 
AGNs with (logt^)L > 9.53 the L[Olll] is almost constant, as 
indicated by its median values in all galaxy luminosity bins. AGN 
hosts with younger stellar populations have the [O III] luminosity 
decreasing with age (see also Kauffmann et al. 2003c, for similar 
results). In other words. Fig. |9l shows that even AGN hosts have a 
bimodal distribution with old galaxies having [O III] luminosities 
lower and constant along age bins (but increasing with galaxy lu- 
minosity), and galaxies with young stellar populations having high 



values of L[0 III], which increases as a galaxy becomes younger 
and brighter. 

In order to investigate this apparent bimodal character of AGN 
hosts, we have divided this class in three subclasses according to 
their emission line properties. The approach used here is analo- 
gous to that of Brinchmann et al. (2004). The subclasses are: (i) 
'composite' galaxies, which appear in t he BPT dia gram between 
the cu rves defined by iKauffmann et al.1 i2003d) and lKewlev et alJ 
i200lh : (ii) hosts of 'Seyfert 2' nuclei types, located in the BPT 
diagram above the curve defined by Kewley et al.; and (iii) hosts 
of low-luminosity AGN (LLAGN), which do not appear in this di- 
agram as they have only the [N ll] and Ha emission lines mea- 
sured with S/N > 3, and are c lassified by using only the crite- 
rion log([N ll]/Ha) > -0.2 fsee lMiller et al.l2003h . In Fig.lTolwe 
show the distribution of light- weighted mean stellar ages for these 



10 Mateus et al. 







o 



AGN hosts 

log L/Lq 

10.07 - 10.29 ^ 

10.29 - 10.49 ": 

10.49 - 11.^6:^ 



18586 
-0.22 




7 8 9 

<iog t^>L [yr] 

Figure 9. Luminosity of [O iii]A5007 emission line as a function of the 
mean light- weighted stellar age. The median values of L[0 ill] as a function 
of age for three ranges in galaxy luminosity are shown as lines with different 
shapes. The vertical dashed line is the age value with best distinguish star- 
forming and passive galaxies. 
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Figure 10. Distribution of mean light- weighted stellar age for distinct types 
of AGN hosts (see text for details). 



three subclasses of AGN hosts, as well as for all galaxies in this 
spectral class. It is interesting to note that composite galaxies at 
one side, and LLAGN hosts at another, are in the tails of the distri- 
bution for all AGN hosts. The majority of AGN hosts in our sam- 
ple is composed by galaxies with old stellar populations, mainly 
represented by the LLAGN hosts which show low values for the 
[O lll]A5007 luminosity. We conclude from the L[0 III] versus age 
relation shown in Fig|9|and the distributions seen in Fig. fTOI that 
the AGN activity is closely linked to the star-formation activity of 
a galaxy, in the sense that galaxies with younger stellar populations 
tend to have more powerful AGNs. 



4.2.2 Stellar mass 

The distribution of stellar mass for our sample is presented in 
Fig.[n] As in Fig.jSj we show the distributions for the full sample 
and for each spectral class discussed here. The value of stellar mass 
which best divides star-forming and passive galaxies is shown as a 
vertical dashed line, and corresponds to ~ 4. 7x10^° M© . The 
reliability and completeness parameters obtained with these values 
are 71 per cent. 



This value is susbtanti ally larger than the 3 x 10^° M© at 
which lKauffmann et alJ (|2^3b ) find a sharp transition in the phys- 
ical properties of galaxies. As mentioned in SEAGal I, this differ- 
ence cannot be attributed to the differences in IMF. Kauffmann et 
al. used the Kroupa IMF, whereas we use the one by Chabrier, but 
as show in the BC03 paper, these two IMFs yield identical M/L 
ratios. 

A more relevant source of discrepancy is related to the ap- 
proach used to determine the transition mass. In tKauffmann et all 
i2003bh . there is no objective criteria to define the stellar mass 
value at which galaxy properties change, whereas we have used 
a clear procedure to find the stellar mass transition at ^ 
4.7 X 10^° M©. 

We have also verified if the use of a volume-limited sample 
in such analysis, with an absolute magnitude limit cutoff M(r) — 
—20.5, affects the determination of the transition mass value, being 
responsible for the difference discussed above. In order to investi- 
gate this issue we have determined the transition mass for galaxies 
in a flux-limited sample containing 20000 galaxies selected at ran- 
dom from the main galaxy sample of SDSS. This results in an even 
larger transition mass, ^ 6.3 x 10^° M©. 

As anticipated in section 12.11 the inclusion of very low Z 
SSPs in the base inevit ably leads to larger stellar masses. The 
iKauffmann e tll] J2003d) mass estimates are based on a library of 
model galaxies constructed with Z ^ 0.25Zq, whereas the SSPs 
in our base go down to O.OOSZ© . This is a systematic source of dis- 
crepancy between the mass scales in these two studies. In SEAGal 
I, which used only Z ^ 0.2^© SSPs, we have shown that the stel- 
lar masses estimated via the spectral synthesis approach are about 
0. 1 dex larger than that obtained by Kauffmann et al. (MPA/JHU 
group). Repeating this comparison for our new stellar masses yields 
a median offset of 0.2 dex. As discussed in SEAGal I, part of this 
difference (~ 0.1 dex) could be attributed to technical details em- 
ployed to estimate the stellar masses. The remaining difference is 
due to the inclusion of low-Z SSPs in our spectral base. We have 
confirmed this hypothesis through the analysis of the results ob- 
tained with a base excluding the 0.002 and 0.005 Z© SSPs. As 
expected, the values obtained by this way are only about 0.1 
dex larger than the stellar masses obtained by the MPA/JHU group. 
The transition mass for this base is ^ 3.9 x 10^° M©. 

Therefore, discrepancy b etween the values of the s tellar mass 
transition obtained here and in IKauffmann etallj2003d) is due to a 
combination of differences in the methodology to define this mass 
and differences in the metallicities in reference set of model stellar 
populations. 



4.2.3 Stellar extinction 

Another important quantity related to physical processes occurring 
in galaxies is the amount of dust present in the interstellar medium. 
Our spectral synthesis approach estimates this value by obtaining 
the attenuation of starlight by dust parametrised in V-band, i.e. the 
stellar extinction. Ay. We plot in Fig.[T2lthis parameter against the 
Dn(4000) index. The results for all galaxies in our sample are de- 
picted in Fig.fT2t. There is a clear anti-correlation between Ay and 
Dn(4000) with rs = —0.66. However, the relation between these 
two quantities is better appreciated when we consider separately 
star-forming, passive and AGN host galaxies. Indeed, star-forming 
galaxies fFig.[T2b) do not show significant correlation between Ay 
and Dni^OOO), with rs = 0.03, and passive galaxies (Fig. IT2l) 
exhibit only a small anti-correlation (rs = —0.17). The median 
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Figure 12. V-band stellar extinction versus Dn(4000) for (a) all galaxy sample, (b) star-forming, (c) passive, and (d) AGN host galaxies. The median values 
of Ay in bins of Dn(4000) index for three luminosity intervals are shown as different lines. 
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Figure 11. Distribution of stellar mass for all galaxies in our sample, and 
for classes of star-forming, passive and AGN hosts. The vertical line is the 
optimal separator between star-forming and passive galaxies. 



of Ay for star-forming and passive galaxies are 0.33 and -0.12, 
respectively. 

The relation between Ay and 1^^(4000) is weak for star- 
forming and passive galaxies, but strong for galaxies harbouring 
an AGN, as can be appreciated in Fig.[T2h. The Spearman correla- 
tion coefficient between these two quantities is —0.60. It is not easy 
to interpret this finding, considering that the measured spectrum of 
this kind of object is a composite of the AGN spectrum plus the un- 
derlying stellar population. A better understanding of what is going 
on requires spatially-resolved spectroscopy of the central regions of 
this kind of galaxy. 

It is interesting to verify how the trends above behave if we 
use galaxy colours. This is shown in Fig.[n]for the colour (u — i). 
The trends of stellar extinction with colour are, qualitatively, quite 
similar to those between Ay and Dn(4000), when all galaxies 
are considered. However, in the case of star-forming galaxies, the 
correlation become stronger with colour (the Spearman coefficient 
being 0.33). This is in agreement with the results reported by 
[stasinska et al. (2004) for the nebular extinction of spiral galaxies: 
it tends to increase with galaxy colour. On the other hand, the small 
correlation shown in Fig.[T2l for passive galaxies disappear in the 
plot of Ay versus {u — i) colour. We d iscuss other prope rties of 
extinction in another paper of this series iSodre et all200d) . 



4.2.4 A note about negative extinction 

As discussed in SEAGal I, for both statistical and physical reasons, 
we did not constrain the extinction Ay to be positive. Interest- 
ingly, Ay < galaxies are nearly all passive, in agreement with 
iKauffmann et alJ i2003ah . who found negative extinction primarily 
in galaxies with a large Dn(4000). As a whole, the distribution 
of Ay for passive galaxies indicates that they have essentially no 
dust, as expected given that they are dominated by old stellar pop- 
ulations. Still, the fact that this distribution is centered on slightly 
negative values of Ay deserves some explaining, as it cannot be 
attributed purely to random errors, which would produce Ay = 
on average. Although this has no impact on the main conclusions 
of this paper, we find it useful to open a parenthesis to look at this 
issue more closely. 

We find that one of the factors which is responsibl e for 
this te ndency is a-enhancement. As reported by Sodre et all 
i2005h . our STARLIGHT fits for passive galaxies underpre- 
dict the strength of a-element features, which is not sur- 
prising given that the STELIB -based BC03 SSPs do not ac- 
count for the a-enhancement known to occur in massive ellipti- 
cals ( Worthev, Faber, & Gonza lez 1992; Davies, Sadler, & Peletiej 
11993: IThomas. Maraston. & BendeJ2002t) . Furthermore, the resid- 
uals in these features are larger as increases, in analogy with the 
[Mg/Fe]-cr^ anti-correlation. 

Gomes (2005) carried out simulations to test the effects of en- 
hancing a-features in test galaxy spectra constructed with the BC03 
SSPs. To emulate this effect realistically, he has added to the mod- 
els the mean residual spectrum of passive galaxies in the regions of 
CN, Mg and NaD. These controlled experiments showed that in or- 
der to try to match these artificially enhanced bands, STARLIGHT 
overpredicts the strength of the older components, where these fea- 
tures are deeper. Furthermore, even though the test galaxies were 
constructed with Ay — 0, the fits yield systematically negative 
extinction, of order Ay ~ —0.1, as found in passive galaxies 
(Fig. [T2b). This happens because the older components invoked 
to match the deeper a- spectral features are also redder, implying a 
continuum mismatch which is compensated by a slightly negative 

A*y. 

The bottom line is that negative extinctions in passive galaxies 
are at least partly due to the use of a base of SSPs which do not 
have the proper chemical mixture. Recently publis hed a-enhanced 
libraries should hel p fixing this mismatch fe.g. ^Coelho et al 120051: 
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Figure 13. Same as Fig.[T2| but now with {u — i) colour in the ordinate. 



5 DISCUSSION 

5.1 Origin of the bimodal distribution 

An easy way to visualise the bimodality of the galaxy population is 
in the colour-magnitude diagram, which shows two separate colour 
seq uences, a broad bl ue sequence, and a narrower red one (see 
e.g. lBaldrv et alfcOO^h . Here we investigate the nature of bimodal- 
ity through the mass-luminosity relation. Fig. 1141 shows the stellar 
mass as a function of the dust-corrected z-band luminosity (in so- 
lar units) for the full galaxy sample and for the galaxy populations 
divided by using the value Dn(4000) = 1.67. Robust linear fits to 
the mass-luminosity relation for these two galaxy sequences yield: 



log(M./M0) 
and 
log(M./M0) 



-4.89 + 1.48 log(L;,/L0), for late-types 



-1.50 + 1.18 log(L^/L0), for early-types. 



Here, the two galaxy sequences commonly seen in the colour- 
magnitude diagram are well characterised with that corresponding 
to early-type galaxies being less steep than that of late-type galax- 
ies. The relation for early-types also shown the best correlation co- 
efficient {rs — 0.93) and less scatter. For late-type galaxies, the 
scattering observed in the — Lz relation is associated to the 
large spread in the mass-to-light ratio and star formation histories 
shown by thes e objects. 

Recently, B aldrv etall i2004h have found that a change in the 
distribution of galaxy colours occurs around (1 - 3) x 10^° M0, 
close to the value at wh ich the galaxy properties also change (see 
iKaufFmann eTal]l2003d) . However, from Fig.[l4|we note that an- 
other parameter is needed in order to achieve a better separation 
between the two galaxy populations characterised by the extreme 
classes of star-forming and passive galaxies. For instance, the value 
log M^/Mq > 10.5 can select almost all early-type galaxies, but 
a large number of massive late-type galaxies is also selected by us- 
ing this criterion. In fact, as shown in Table |5| the reliability and 
completeness parameters for star-forming and passive galaxies ob- 
tained with the optimal value of log M^/Mq = 10.67 are lower 
than those obtained when considering the mean light- weighted stel- 
lar age to split the galaxy populations. Thus, it seems that the galaxy 
mass is not playing a main role in defining the bimodal galaxy dis- 
tribution seen in local galaxies, reinforcing the idea that the ob- 
served bimodality is related to the presence of a young component 
in low-mass, blue galaxies. 

Another way to investigate this issue is presented in Fig. fTsl 
where we show the {u — r) colour and Dn (4000) index versus the 
mean light-weighted stellar age and stellar mass for star-forming 
and passive galaxies in our sample. The values which better divide 
these extreme galaxy populations are shown as vertical and hor- 



izontal dashed lines. The median values of colour and Dn(4000) 
for three bins of galaxy luminosity are also shown as different lines. 

With the help of this figure, we have investigated the trend 
of the mean light-weighted stellar age to distinguish the extreme 
galaxy populations in a better way than if using stellar mass. Con- 
sidering the optimal values of {u — r) colour and Dn(4000) index, 
and that of (logt^)^ and logM^, to characterise early and late- 
type galaxies, there are two sets of abnormal galaxies (located in 
the top-left and bottom-right quadrants of Fig.lTsl: late-type galax- 
ies (e.g. with Dn(4000) < 1.67) with old stellar ages or high stel- 
lar masses, and early-type galaxies (with Dn(4000) > 1.67) with 
young stellar ages or low stellar masses. We find that the fraction 
of late-type 'old' plus that of early-type 'young' galaxies is about 
6 per cent, whereas the fraction of late- type 'massive' plus early- 
type 'low-mass' galaxies is 29 per cent. Thus, taking a single value 
to characterise the two main galaxy populations (star-forming and 
passive), for instance based on the Dn(4000) index, the fraction 
of those uncommon galaxies is higher when we consider the opti- 
mal value of the stellar mass, in contrast with the fraction obtained 
when considering the mean stellar age to split the galaxy popula- 
tion. A similar trend is also seen when all galaxies are taken into 
account. 

Thus, supported by these findings, we argue that the bimodal- 
ity of the galaxy population, commonly seen in colour-magnitude 
diagrams, single colour distributions, and in the mass-luminosity 
relation discussed above, is related to the presence of a young and 
luminous stellar component in galaxies currently undergoing star 
formation, in contrast with the stellar content of passive galaxies, 
where old stars are the responsible for most of their luminosities. 



5.2 The role of stellar mass 

Colours of galaxies reflect their star formation histories, which are 
also linked to the stellar populations that we are seeing today. For 
instance, a galaxy that has suffered recent bursts of star formation 
will contain a considerable fraction of young and hot stars that will 
be responsible for a large fraction of the light we receive from it. 
On the other hand, earlier events of huge star formation activity will 
keep the fraction of old stellar populations even higher, leading to 
galaxies with redder colours and spectra dominated by absorption 
lines. 

We also observe that red, passive galaxies tend to be more 
massive than blue, star-forming galaxies. Thus, the stellar mass of 
a galaxy is related to the bimodality seen in the colour distribution 
of local galaxies, in the sense that only less massive galaxies appear 
to be forming stars in the last few gigayears. 
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Figure 14. Mass-luminosity relation for (a) all galaxies in our sample, and for (b) late-type galaxies with Dn(4000) < 1.67 and (c) early-type galaxies with 
Dn(4000) > 1.67. Solid lines are robust fits for the relations. 
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(logt^>L 


9.53 


94.8 


92.9 
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92.0 


77.3 


log M^/Mq 


10.67 


64.2 


79.9 


83.0 


59.2 


25.2 



Table 5. Reliability and completeness parameters for concentration index, (u — r) and (u — i) colours, Dn(4000) index, mean light- weighted stellar age and 
stellar mass optimal values obtained to split the distributions of star-forming and passive galaxies. 



We have investigated this tendency by examining the mass- 
weighted mean stellar age, (t^)M ^, of galaxies in our sample, as 
a function of their stellar masses. This age is related to the epoch 
of formation of the stellar population which nowadays contributes 
significantly to the galaxy mass. In Fig.[l6lwe show the relation be- 
tween stellar mass and (t^)M (in Gyr), for (a) all galaxy sample, (b) 
star-forming, and (c) passive galaxies. There is a clear correlation 
between stellar mass and (t^)M, with a Spearman rank coefficient 
of rs = 0.69 for all galaxy sample. The median values of stellar 
mass in bins of stellar age containing the same number of objects 
are also shown in this figure. Note that the median stellar mass in- 
crease from 1.4 X 10^° to 1.5 x 10^^ M©, from the young to the 
oldest age bin. This trend is particularly related to the existence of 
a dominant fraction of passive galaxies, with high stellar masses 
and without signs of recent star formation activity, at the oldest age 
bins shown in Fig.[T6t. In addition, less massive galaxies are still 
forming stars which will contribute to increase their masses further. 

When examining these plots, it is fit to recall that we are us- 
ing a volume-limited sample. The limiting magnitude of M(r) = 
—20.5 (corresponding to M* (r) + 1) is large enough to detect mas- 
sive galaxies of any age, but old low-mass systems will be excluded 
by construction. To examine this issue quantitatively, we use the 
SSPs in BC03 and determine, for each age and metallicity, the stel- 
lar mass which results in M(r) = —20.5. The results are shown 
as solid lines, one for each Z, in Fig. [n] overplotted to the 
versus (t^)^ and versus {t^)M relations for our sample. The 
M^(t, Z)M(r)=-20.5 lines trace very well the lower envelope of 
the data points in the left panel of Fig. [n] The lower envelop be- 
comes fuzzier using (t^)M, but the overall conclusion is still the 



1 actually {U)m = lO^^^g**)^ 



same. Except for a handful of cases, all galaxies in this sample lie 
above these lines. Old, low mass objects ("dwarf elliptical"-like) 
which would populate the lower-right region of this diagram are 
too faint to satisfy our M(r) < —20.5 cut. On the other hand, and 
perhaps more interestingly, the nonexistence of massive galaxies 
with young stellar populations (mainly noted in the relation with 
(t*)M), which should appear at the top-left corner of Fig. [n] is 
clearly unrelated to our sample selection. Thus, our claim that mas- 
sive galaxies have essentially older stellar populations whereas, in 
contrast, less massive ones have younger populations, does not de- 
pend on the limiting magnitude adopted in this work. 

5.3 Downsizing in galaxy formation 

In a biased galaxy formation scenario fe.g. lCen & OstrikeJllQQBh . 
massive galaxies form from the highest peaks in the initial den- 
sity fluctuation field. In this scenario, massive galaxies formed at 
high redshifts tend to inhabit the high-density region s associated to 
galaxy clusters and rich groups seen today. Indeed, Ko dama et alJ 
i2004) in a photometric analysis of galaxies in colour- selected 
high-density regions at 2; ^ 1, have shown that galaxy forma- 
tion processes, including mass assembly and star formation, take 
place rapidly and are completed early in massive systems, while in 
the less massive ones these processes (at least star formation) are 
slower. 

This galaxy formation scenario, in which more massive 
galaxies are formed at higher redshift, was first proposed by 
.Cowie et al.. Il996ii) . who found that the maximum luminosity (or 
mass) of galaxies undergoing rapid star formation has been declin- 
ing smoothly with decreasing redshift since z ^ 1. Recently, many 
observational results have supported this idea by suggesting that the 
most luminous galaxies formed the bulk of their stars in the first 
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Figure 15. Mean stellar age and stellar mass as a function of (u — r) colour (top panels) and Dn (4000) index (bottom panels), for star forming and passive 
galaxies in our sample. It is also shown the number of objects and the Spearman rank correlation coefficient for each panel. 



^ 3 Gyr of cosmic hist ory iMcCarthv et alJl2004l : iKodama et alJ 
|200i;|Juneauetal|2005|, and others). 

Hence, even z ^ 1 there would be a red and luminous 
population composed by massive and old galaxies, and a blue 
and faint one, formed by less massive galaxies with still ongoing 
star formation. Recently, evidence for this bimodal distribution of 
galaxy properties (at least in colour) around this redshift has been 
found by many works throug h different approaches and using dif- 
ferent observational data sets JBell et all2004l : IWiegert et all2004l : 
iWeiner et al. 2005). 

In the local universe, the bimodality of galaxy properties 
was firstly noted in the colour distribution of SDSS galaxies 
(e.g. IStrateva et alJl200lT) and in the star formation properties of 
galaxies from the 2dFGRS (e.g. Madswick et al. 2002). In the 
last years many works have focused their analysis on the de- 
pendence of galaxy prop erties on the stellar mass. For instance, 
iKauffmann et alJ i2003d) have shown that at a stellar mass above 
3 X 10^° M© there is a rapidly increasing fraction of galaxies with 
old stellar populations, high surface mass densities and high con- 
centrations typical of bulges. These results, and similar ones based 
on th e star formation properties of galaxies (e.g. lBrinchmann et alJ 
l2004h . together with our own findings derived from the analysis 
of the mass-weighted mean stellar age of galaxies, are in perfect 



agreem ent with the 'downsizing' scenario proposed bv lCowie et alJ 

Notwithstanding the numerous observational reasons to be- 
lieve in this scenario, it is still very puzzling when we assume 
the hierarchical galaxy formation (bottom-up) scenario of ACDM 
models, which predicts that small systems collapse first and mas- 
sive galaxies form later via the assembly of these already col- 
lapsed structures. The theoretical 'downsizing' picture still re- 
quires significant improvements, mainly concerning galaxy forma- 
tion and related proc esses occurred since z = 1, as pointed out by 
iHammer et alJ I i2005h . At this point it is also important to stress the 
recent advance of semi-analytic models based on hierarchical clus- 
tering to explain some observational properties of galaxies, such 
as the bimodality in galaxy colours (e.g. Menci et al . 2005). In this 
framework, distinct merging and star formation histories for pro- 
genitors with different masses lead to a high- 2; origin of the bi- 
modality seen in the properties of galaxies up to z 1, particularly 
during the formation and merging of their progenitors. 

Despite the uncertainties in the scenarios for galaxy evolution 
and the current status of 'downsizing' in galaxy formation and hi- 
erarchical models, here we have shown that our spectral synthesis 
method discussed in SEAGal I enables us to investigate the star for- 
mation history of local galaxies, making a link among these objects 
and their progenitors, which have been recently explored by cur- 
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Figure 16. Stellar mass versus the mass-weighted mean stellar age of formation of the stellar population which more contributes in mass, for (a) all galaxy 
sample, (b) star-forming, and (c) passive galaxies. The solid lines are the median and quartile values of stellar mass in age bins containing the same number of 
galaxies. 




5 10 15 5 10 



<t.>L [Gyr] <t,>^ [Gyr] 

Figure 17. Stellar mass as a function of the mean Hght-weighted stellar age (left panel) and mass-weighted stellar age (right panel). The soHd lines represent 
the stellar mass for which a SSP of age (t^) l has an absolute r-band magnitude of —20.5. Each line corresponds to a stellar metallicity. From bottom to top: 
Z = 0.005, 0.02, 0.2, 0.4, 1 and 2.5 Zq. 



rent high-redshift galaxy surveys. In future works we will discuss 
in detail this and other related issues. 



6 SUMMARY 

In this second paper of the SEAGal collaboration, we have investi- 
gated the bimodality observed in galaxy populations by inspecting 
the spectral properties of SDSS galaxies. We have used the physi- 
cal parameters derived from the spectral synthesis method applied 
to a sample of about 50 thousand galaxies extracted from the SDSS 
Data Release 2. Galaxies are classified according to their emission 
line properties in three distinct groups: star-forming, passive, and 
AGN hosts. The bimodality of galaxy properties is investigated 
with emphasis on these spectral classes. Our main findings are sum- 
marised below: 

(i) The bimodality of the galaxy population can be represented 
by two extreme spectral classes, corresponding to star-forming 
galaxies at one side, containing young stellar populations and pref- 
erentially of low stellar masses, and passive galaxies at the other 
side, without ongoing star formation and populated by older stars. 

(ii) In an intermediate locus, there are the AGN hosts, which 



comprise a mix of young and old stellar populations. However, this 
spectral class also shows a bimodal behaviour when considering the 
[O III] emission line luminosity, with younger galaxies showing the 
larger values of L([0 III]), and older ones presenting the same low 
luminosities {L{[0 III]) < 10^ L©) essentially at all stellar ages. 

(iii) As a main result, we found that the mean light-weighted 
stellar age of galaxies is the direct responsible for the bimodality 
observed in the galaxy population. 

(iv) The stellar mass, in this view, has an additional role since 
most of the star-forming galaxies present in the local universe are 
low-mass galaxies. Our results also reinforce the idea of a 'down- 
sizing' in galaxy formation, where massive galaxies seen nowadays 
have stopped to form stars more than 10 Gyr ago. 



In this work we have explored a piece of the arsenal of physi- 
cal parameters provided by the spectral synthesis method, in order 
to investigate the spectral properties of galaxies and revisited the 
bimodal character of the galaxy population. Other papers in this 
series will bring more insight to some fundamental problems con- 
cerning to the physical properties of local galaxies. 
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APPENDIX A: APERTURE EFFECTS 

A particular effect acting in modern galaxy redshift surveys is 
related to the fixed- size apertures used to obtain spectroscopic 
data. In the 2dFGRS the fibre diameter is only 2 arcsec, whereas 
in the SDSS the fibres cover 3 arcsec of the sky, leading to a 
small sampling of the integrated light of nearby galaxies. Recently, 
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Figure Al. (Top) Cumulative distribution of the ratio of the fibre radius 
to the galaxy radius for all sample and for each spectral class discussed 
in this work. (Bottom) Cumulative distribution of the correction factor to 
be applied to correct for aperture effect, also for all sample and spectral 
classes. The vertical dashed lines in these plots are the median values of the 
distributions for all sample. 



iKewlev. Jansen & Gelleil J2005h have examined this problem in de- 
tail by investigating the effect of aperture size on the star formation 
rate, metallicity, and extinction determinations for galaxies selected 
from the Nearby Field Galaxy Survey. The main result obtained 
by these authors is the need to select galaxies with z > 0.04, in 
the case of the SDSS, to minimise aperture effects in the spectral 
measurements. This redshift limit is required to assure that a fibre 
captures more than 20 per cent of the galaxy light. 

In this work, the galaxy sample was built to contain only 
galaxies with redshifts larger than z — 0.05. Hence, following the 
recommendation given by Kewley et al., the aperture effects in our 
sample tend to be strongly minimised. On the other hand, in Fig.0 
we have shown that the fraction of star-forming galaxies increases 
significantly with increasing redshift, probably due to the fact that 
in nearby galaxies the fibre is capturing mainly their bulge regions, 
whereas the nebular emission lines (used here to classify galaxies) 
come essentially from their discs. Here, we inspect some aspects of 
this effect. 

In Fig. lAll we show the cumulative distributions of two quanti- 
ties related to the covering fraction of the fixed- size aperture for all 
galaxy sample, and for each spectral class discussed in this work. 



In the top panel of this figure, it is shown the ratio between the fi- 
bre radius (1.5 arcsec) and the Petrosian radius containing 90 per 
cent of the galaxy flux (in r-band). Thus, this ratio represents the 
fraction of the galaxy size that is actually covered by the fibre. The 
median value for all galaxy sample, and also for each spectral class 
taken individually, is about 23 per cent. Thus, the fibre covers, on 
average, the same fraction of galaxy size almost independently of 
the galaxy type. 

The bottom panel of Fig. lATl shows the aperture correction. A, 
to be applied in flux related measurements, based on the differences 
between the total galaxy magnitude in the r-band and the mag- 
nitude inside the fib re (r fibre), explicitly A = iQ-^-'^i^-'^ fibre) 
taopkins et aljEoOSi) . This correction factor is commonly applied 
to emission line derived parameters (like Ha SFR), and assumes 
that the regions where the lines are formed (mainly H II regions) 
have a distribution identical to that where the continuum is pro- 
duced by the old stellar populations. In Fig. lAll we note that 
star-forming galaxies have aperture corrections larger than that of 
other classes. This may be somewhat problematic, because this 
kind of correction i s performed just in this type of galaxies (e.g. 
iHookins et al.l2003h . 

From the plots shown in Fig. lAll we conclude that while the 
3 arcsec diameter fibres cover about the same fraction of galaxy 
size for all galaxies, the fraction of light gathered by their small 
apertures are biased towards lower values (or consequently higher 
aperture corrections) for star-forming galaxies, and higher values 
for passive galaxies. Since the expected aperture effects would tend 
to increase the number of passive galaxies relative to that of star- 
forming at lower redshifts, the low aperture corrections obtained 
for these galaxies indicate that, on average, there are not substantial 
effects acting on the fraction of this spectral class along the redshift 
range considered in this work. 

Al Mass- to-light ratio 

In order to estimate the total stellar mass of a galaxy from the spec- 
tral synthesis, we have assumed that it has a constant radial M/L. 
In other words, total stellar masses are computed from the fibre 
and Petrosian (total) magnitudes by considering the M/L inside 
the fibre as representative of the whole galaxy. The effects of this 
assumption on our stellar mass determination is investigated here. 

We have carried out the following tests to address this inter- 
esting question. From the SDSS photometry data we have both fi- 
bre and total fluxes and colours, from which one can also compute 
outside-the-fibre fluxes and colours. We have examined the corre- 
lations between inside-the-fibre M^/Lz and fibre-colours and cho- 
sen the M^/Lz versus (r — z) colour as the best one (though other 
colours yield about equally good correlations). Having established 
an empirical relation parametrised by M^/Lz = a{r — z) -\- b 
(which we did separately for passive and star-forming galaxies), 
we can then, from an observed outside-the-fibre colour estimate 
the corresponding outside-the-fibre M^/Lz. 

In the top panel of Fig. IA2I we show the relation between 
the M^/Lz obtained inside the fibre and outside it, for star- 
forming and passive galaxies. Compared to the inside-the-fibre 
value, M^/Lz is, in the median 0.13 dex smaller in the outside. 
This applies to star-forming galaxies. This is expected, as the out- 
side light is dominated by disk-light, whereas the fibre data also 
includes the bulge. Differences for the passive galaxies were negli- 
gible (-0.03 dex). 

Total masses computed with this approximate correction are, 
in the median, only 0.09 dex smaller than those obtained extrapo- 
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Figure A2. Top: relation between the mass-to-light ratio obtained inside the 
?)" fibre and outside it (see details in text) for star-forming galaxies (blue 
points) and passive galaxies (red points). Bottom: distribution of the differ- 
ences between the total stellar masses estimated by considering a constant 
/Lz along the radius of a galaxy and estimated via the outside-the-fibre 
M^jLz, for star-forming galaxies (blue) and passive ones (red). 



lating the M^/Lz derived from the fibre- spectroscopic data to the 
whole galaxy, as can be seen in the bottom panel of Fig. lA2l Be- 
cause this is a small difference, and also because the M^jLz versus 
colour calibration has a good deal of scatter, we prefer not to apply 
this correction in our analysis. 
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author. 



